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ABSTRACT 


With continuing improvements in rare earth permanent magnet (REPM) technol- 
Ogv, applications for their use are being discovered that were previously not possible. 
Two such applications for permanent magnets are in focusing elements for linear accel- 
erators and ion sources, and in insertion devices (wigglers and undulators) used to 
produce synchrotron radiation. However, these magnetic transport elements are sub- 
jected to high radiation levels. Consequently, there is considerable interest in the United 
States and abroad to discover and quantify the effects of radiation on REPMs. Using 
the Lawrence Livermore National Laboratory (LLNL) 100 MeV Linac, four different 
samples of REPM were irradiated to one to two gigarads of exposed dose from a 
bremsstrahlung production target in an attempt to simulate the consequences of beam 
spilis of a high energy primary electron beam. Of the samples irradiated, Sm,Co,, proved 
to be the most resistant to gamma radiation. The electron transport code CYLTRAN 
of the Integrated Tiger Series (ITS), which is an electron and photon Monte Carlo sim- 
ulation code, was used to determine the angular and energy spectra for both electrons 


and photons produced by the target used at the LLNL Luinac. 
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I. INTRODUCTION 


The development of permanent magnet materials has progressed tremendously since 
the beginning of this century. However, new to technology then, permanent magnets 
gave a relatively low performances in most technological applications. As a result, 
electromagnets were preferred. With steady increase in quality and performance, per- 
manent magnets are supplanting electromagnets for many applications. Today perma- 
nent magnets are so widely used that one could consider them as an indispensable 
prerequisite of modern technology. By far the largest application of permanent magnets 
is In motors and generators. Next in importance are telecommunication, data manipu- 
lations and for measurement and control applications. Substantial numbers of perma- 
nent magnets are also used in acoustic devices and magnetomechanical applications. 
cet. I: p. 11] 

The application depends on the quality or strength of the permanent magnet. A 
measure of quality for permanent magnets is the maximum energy product (BH), 
occuring in the magnetic hysteresis loop. The maximum energy product is commonly 
expressed in megagauss-oersteds (MGOce). In the SI unit system, (BH),,,, 1s measured 
in kJ m°, where one {GOe equals 7.96 kJ’m*. An ideal permanent magnet would have 
a constant magnetization equal to the saturation magnetization M, in anv applied field. 
If pure iron could be made into an ideal permanent magnet with 100 percent density, it 
would have an energy product of 107 MGOe [Ref. 1: p. 11]. This is considered to be the 
maximum energy product achievable in a magnetic material. Appendix A provides de- 
tailed descriptions of magnetic properties of materials and Appendix B is a glossary of 
common terms associated with magnetic properties of materials. 

Since the mid 1960’s, improvements in rare earth permanent magnet (REPM) tech- 
nology have shown remarkable progress [Ref. 2: p. 1]. Figure I on page 2 shows the 
historical development of the energy product quantity. Improvements in the quality of 
compounds of rare earth materials such as samarium-cobalt (Sm-Co) and more recently 
compounds of neodymium, iron and boron (Nd-Fe-B) have provided new applications 
that were previously not possible. The use of REPMs in accelerator applications include 
magnetic focusing devices for linear accelerators and ion sources, and insertion devices 
(wigglers and undulators) for synchrotron radiation sources. These applications have 


evoked considerable interest in permanent magnets both in the United States and 
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Figure 1. Improvements in Magnetic Quality of Permanent Magnets 


abroad. These materials are characterized by a high remanent induction B, and high 
coercive force H, which are used to produce high multipole fields in extremely compact 
magnets. These applications often require the magnets to operate in high radiation en- 
vironnients. 

Although manv magnetic properties of REPM materials have been studied, to date 
the understanding of the effects of ionizing radiation is phenomenological at best. Pre- 
vious studies have focused mostly on determining the effects of exposing REPM samples 
to beams of protons and fission neutrons. However, the permanent magnet components 
of a wiggler in free electron laser (FEL) systems can be expected to receive very poten- 
tially high doses of secondarv electron and photon radiation. High radiation levels can 
be produced from electrons that are lost from a primary beam due to accidental beam 
spillage, poor emittance, or improper front-end matching. The lost electrons interact 
with the beam pipe producing a mixed field of radiation, including bremsstrahlung 
photons, a fission-like spectrum of neutrons and a broad distribution in energy of 
electrons [Ref. 3: p. 2]. Experiments were conducted at the Linear Accelerator (Linac) 
at Lawrence Livermore National Laboratory (LLNL) to determine which commercially 
available permanent magnet alloys are the most insensitive in their magnetic properties 
to high levels of this tvpe of radiation environment. 

The Integrated Tiger Series (ITS), a Monte-Carlo simulation code developed to 
study the manner in which electrons and photons propagate through materials, was used 
to determine electron and photon energy and angular spectra of radiation produced by 
the target at LLNL’s Linac. The ITS Monte-Carlo scheme calculates a particle's energy, 
direction and position at discrete intervals as the particle travels through a material. 
Energy deposition and ionizing levels are computed. Total dose in any material can then 
be calculated. 

The ITS code CYLTRAN simulates the transport of electrons and photons through 
a cylindrical geometry. The beam production target and experimental setup at LLNL’s 
Linac were modeled using CYLTRAN. Results from the simulation codes were com- 


pared with experimental data and are discussed in Chapter IV. 


Il. BACKGROUND 


In our experiment conducted at the 100 MeV Linac Facility at LLNL, studies were 
made to measure the effects on high levels of gamma radiation on rare earth permanent 
magnets. Other experiments dating back to July, 1982 were conducted by other labo- 
ratory groups. These experiments however, were focused on radiation effects on mag- 
nets exposed to high fluences of protons and neutrons. 

In 1982, radiation effects on Sm-Co permanent magnets were studied at Los Alamos 
National Laboratory (LANL). Samples of SmCo, (Hicorex 90B) and (SmPr)Co, 
(Hicorex 96B) manufactured bv Hitachi Magnetics were irradiated by neutrons produced 
at the Los Alamos Meson Production Facility (LAMPF). The neutrons were produced 
by the interaction of the LAMPF 800 MeV proton beam with various isotope pro- 
duction targets. A maximum neutron fluence of 1.1 x 10'* n/cm? was deposited on each 
magnet. A two percent drop in magnetic remanence was measured after the irradiation. 
Temperatures for the SmCo, and (SmPr)Co, samples were monitored with 
thermocouples positioned near each sample. The maximum temperature recorded was 
140 °C at maximum proton beam current. Since the stabilizing temperature on these 
magnets Were designed at 100 °C, the studies concluded that the decrease in remanence 
Was consistent with thermally induced demagnetization. No clear conclusions were 
made in this experiment to what effects the radiation had on the Sm-Co magnets. [Ref. 
2: pp. 4-6] 

In 1983, studies on proton radiation effects on rare earth cobalt (REC) permanent 
magnets were conducted at CERN. These studies were designed to examine radiation 
effects on the relatively high remanent induction B, and high coercive force H, magnets 
used for particle beam transport elements such as quadrupoles and sextupoles. However 
little was known about the radiation resistance of RECs. Numerous applications of 
REC magnets at CERN made necessary studies in this area. Four different samples of 
Sm-Co were irradiated. The magnets used were two samples of SmCo, (Vacomax 200 
and Koermax 160), one sample of Sm,Co,. and the fourth sample was a composite ma- 
terial composed of 65.6 percent copper, 24.1 percent samarium and 10.3 percent com- 
posite of other various rare earth metals (Recoma 20). The energy of the incident proton 
beam used to irradiate the magnet samples was 400 GeV. A total dose of 11.4 gigarads 


was deposited on the SmCo,; (Koermax 160) sample, and 10.5 gigarads on the Sm,Co,, 


sample. An unexpectedly high loss in remanence was found on all samples with the ex- 
ception of Sm,Co,,.. The Sm,Co,. sample, which received a dose of 10.5 gigarads, suf- 
fered a three percent loss in remanence. The SmCo, (Koermax 160) sample received a 
total dose of 11.4 gigarads and suffered a much greater loss in remanence, 24.2 percent 
[Ref. 4: p. 10]. Thermal effect studies were performed prior to irradiation to demonstrate 
that a loss in remanence was not caused by excessive temperature. Magnet samples were 
heated to 90 °C for 120 hours in an oven and the measured magnetization loss was less 
than 0.25 percent. During irradiation, the temperature fluctuated between 30 °C and 
80 °C. According to manufacturer specifications, the demagnetization effect at these 
temperatures was less than one percent. Therefore, it was concluded that there was no 
significant loss due to temperature [Ref. 4: p. 3]. 

A similar study was performed at TRIUMPF, in Vancouver, B.C. in 1985. Like the 
study performed at CERN, Sm-Co permanent magnets were being used extensively in 
radiation environments. Four samples of different Sm-Co compounds and one sample 
of Nd-Fe-B were irradiated with a 500 MeV proton beam. Hitachi Magnetics Corpo- 
ration supplied samples of SmCo, (Hicorex 90B) and (SmPr)Co, (Hicorex 96B). Crucible 
Magnetics supplied another sample of SmCo, (Crucore 18) and Sm,Co,- (Crucore 26). 
The Nd-Fe-B samples were supplied by IG Technologies (Ne[GT 27). After receiving 
a total proton dose of 12 gigarads, Crucore 26, a Sm,Co,, sample proved to be the most 
resistant to radiation, losing less than one percent of its original magnetic remanence. 
Sm.Co,- was found to be 10 times less sensitive to proton radiation compared to 
SmCo,. consistent with data measured at CERN. The Nd-Fe-B samples proved to be 
extremely sensitive to ionizing radiation, losing essentially all of its magnetization at a 
dose of 70 megarads and over 50 percent loss at four megarads, more than five order of 
magnitudes more sensitive than Sm,Co,; [Ref. 5: pp. 1-5]. Thermal effect studies were 
conducted using the same techniques utilized at CERN, except that the samples were 
heated to a temperature of 250 °C. The Hicorex 90B suffered the greatest magnetization 
loss from thermal effects, a two percent drop in remanence. The maximum temperature 
reached during radiation was 125 °C for a period of 30-75 hours. Therefore, temperature 
was not considered the major factor in remanence loss (Ref. 5: pp. 1-4]. 

In 1986, LANL conducted a further study of radiation effects on Nd-Fe-B and 
Sm-Co magnets at the Omega West Reactor Site. Gamma radiation effects were also 
studied on two different samples of Nd-Fe-B using a Co® source. However, the Nd-Fe-B 


samples that were irradiated by gammas to a total dose of 49 megarads exhibited no 


decrease in remanence loss. [Ref. 6: p. 3]. This result coupled with the proton 
irradiation measurement from TRIUMPF, suggest that for Nd-Fe-B, the loss of 
remanence is not strictly a total dose effect, but radiation mechanism dependent as well. 
Nd-Fe-B sample proved much more sensitive than the $m,Co,,, which is consistent with 
previous experiments. Also, three different manufactured Nd-Fe-B samples were 
radiated at LANL reactor site with equivalent maximum neutron fluence. Each sample 
measured a different amount of remanence loss, which suggest that radiation effects 1s 
also dependent on the microstructure and manufacturing techniques of the magnetic al- 
lov. Table 1 on page 8 summarizes the results from the previous experiments just 
mentioned. 

The radiation runs scheduled at LLNL’s Linac were directed primarily to continue 
research in gamma radiation effects on Sm-Co and Nd-Fe-B permanent magnets. The 
prior studies established the effects neutron and proton radiation have on these magnets. 
Aside from Co® gamma radiation conducted at Los Alamos, experimental results were 
essentially non-existent. Therefore, a team at LLNL, headed by Dr. Nicholas Colella 
and Dr. Bernhard Kulke, scheduled radiation runs to continue an in-depth study to ex- 
amine the responses of permanent magnets used in FELs due to radiation damage ex- 


posure and consequenct damage. The following 1s a list of magnetic characteristics [Ref 
opal 


]. Halbach hybrid wigglers require REPMs to operate near their coercivity points and 
to possess nearlv linear B-H curves in proximity to ther coercivity points. Non- 
linearities introduced into the B-H curve by radiation damage can profoundly de- 


grade the performance of a wiggler. 


2. A high degree of cancellation of field errors demands that magnets be tightly 
matched in remanence. Further dispersion of remanence caused by radiation 
damage will degrade the initial matching. Hence, variances in remanence of a col- 


lection of magnets must be determined based on post irradiation results. 


3. The radiation fields emploved in previous studies are not representative of the 
mixed fields possible in high powered FELs. Electrons lost from a primary beam 
due to accidental beam spillage, poor emittance, or improper matching, will interact 
with the beam pipe to produce bremsstrahlung photons, a fission-like spectrum 
(sub MeV to several MeV) of neutrons and a broad distribution in energy of 


electrons (deltas and multiple scattered primaries). Thus, in addition to 


6 


nucrostructural damage caused by scattered neutrons and ionization caused by 
neutrons and electrons, REPMs in high powered FELs can suffer total dose 


damage from bremsstrahlung photons as well. 


4. The creation of an alloy with high coercivity is often described as an “Edisonian” 
process [Ref. 7: p. 6]. The elemental composition, particle size and details of 
manufacturing a standard REPM, Sm,Co,, for example, vary widely in industry. 
Thus, one must exercise caution when using the results from one commercial alloy 


to predict the behavior of another in an FEL application. 


The studies at LLNL are considered extensions of the experiments conducted at LANL, 
CERN and TRIUMPF and were performed to provide input to the selection of magnetic 
allovs for future high-powered FELs. The following chapter gives a detailed account of 


the irradiation runs at LLNL’s 100 MeV Luinac. 


Table 1. SUMMARY OF PRIOR STUDIES 


Type of Maximum Dose Remanence 
Sample Alloy Irradiation (x 10° rads) loss (%) 
CERN, 1983: [Ref. 4] 
RECOMA 20 RECO: proton 9.700 42,70 
VACOMAX 200 SmCo, proton 10.400 106.10 
KOERMAX 160 SIO proton 11.400 24.20 
KRUPP WIDIA Sie or proton 10.500 2.60 
TRIUMPF, 1985: [Ref. 5] 
HICOREX 90B Sim or proton 3:02 [3250 
MiCOna sos (SmPr)Co, proton | sie: 6.50 
CRECORE SiGe. proton Sys 1.64 
CRUCORE 26 Si Co. proton 5.94 0.30 
NCIGH 27 cee ab proton 0.003 55.40 
LANL, 1986 (Cobalt 60 Source): [Ref. 6] 
CRU MAN 282 Sa Bae: ganima 48.8 MRads 0.00 
mel Crh 27 INGE ns gamma 48.8 MRads 0.00 
Maximum Fluence 
(x 10 n/cem?) 
LANL, 1982 (LAMPF): [Ref 2] 
HICOREX 90B Sin Oe neutron 1.10 1.88 
HICOREX 96B (SmPr)Co, neutron 1.20 PRIS | 
LANL, 1986 (Omega West Reactor Site): [Ref. 6] 
CRUMAX 282 Nd) Fe).b neutron 2.50 79.10 
NelG i 27m NG Pea neutron 2250) $6.80 
HICOREX 94B Nae ws neutron 3.80 14.00 
INCOR 18 SmGo- neutron 2.60 0.00 
INCOR Zitz $m,Co,, neutron 2.60 0.20 


ih EXPERIMENT 


AAs mentioned in the previous section, studies on radiation effects on permanent 
magnets have focused primarily on effects due to proton and neutron irradiations. Since 
only few experiments existed concerning gamma radiation effects on REPMs, irradiation 


runs were scheduled to expose REPMs to high dosage levels of gamma radiation. 


A. PRELIMINARY CALIBRATION MEASUREMENTS 

Several gamma production targets were available at LLNL’s Linac. The target 
chosen was a composite material known as “heavimet.” This material is composed of 
four percent copper, seven percent nickel and 89 percent tungsten and has a density of 
17.1 gr cm?. Since tungsten 1s the dominant material in this foil, the production target 
will be occasionally referred to as a tungsten target. Figure 2 on page 10 is a cross 
sectional view of the target. The production target is composed of the tungsten foil 
which is attached to a copper layer that forms the lining of the water-cooled support 
Structure. Hitting the production target with an 85 MeV electron beam provided a 
source of ganima radiation. The irradiation runs were conducted in January and Feb- 
ruarv, 1988. 

Prior to the scheduled irradiation runs, preliminary bremsstrahlung dose measure- 
ments were made to determine the radiation dosage pattern from the beam center to a 
radius of 10 cm. This geometry was chosen to accomodate the magnet sample-holder 
used in the experiment. LiF TLD chips were positioned 38.1 cm from the production 
target to measure the dose. The beam energy used for the calibration run was 79 MeV. 
The measured radiation dosage pattern, shown in Figure 3 on page 11[Ref. 8: p. 1], was 
scaled to $5 MeV and 78.7 cm to correspond to the actual experimental parameters. 
Geometrically, the resultant pattern width was scaled as R and the dose was scaled as 
1/R?. The measured on-axis dose was scaled as (MeV)**, a conversion scaling factor 
derived from a flash x-ray calculation. An on-axis dose rate of 0.201 gigarads per 
coulomb of primary electron beam charge was calculated. [Ref. 9: p. 3] 

One of the primary concerns of the Linac operator was ensuring that the beam was 
focused to a size small enough to ensure complete capture by the tungsten disc. To de- 
tect whether the beam diameter was within the disc, glass plates were exposed at the 
TLD chip location to obtain a continuous intensity profile. If a halo effect (beam di- 


ameter greater than the tungsten disc) was present, the glass darkened due to its 


Heavimet 





Figure 2. Longitudinal Cross Section of the Production Target: where the 
heavimet disc has radius of 0.94 cm and the copper support structure has 
a radius of 3.94 cm. 


exposure to beam electrons. The radiation dosage pattern in Figure 3 was calibrated 
assuming that no halo effect existed. The beam radius was estimated at 0.75 cm. [Ref. 
8: pe )] 

Prior to the irradiation, Lif TLD chips were used to verify the gamma dose cali- 
bration measurements taken in April, 1987. An array of TLDs were positioned 167.8 
cm from the production target. The dosimeters were placed at the corners and behind 
a 10 cm lucite block. The center of the bloc!. was positioned on the beam axis. A total 
charge of 11 microcoulombs were deposited. The dosimeters measured 22 +5 megarads 


per coulomb. The pre-calibrated dosage pattern should have given a rescaled dose of 
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Figure 3. Normalized Radiation Pattern of Production Target 


37 =15 megarad per coulomb. The disparity could have been generated from the con- 
version factor derived from the flash x-ray scaling. Keep in mind the large uncertainties 
produced from the flash x-ray conversion factor [Ref. 9: p. 8]. These parameters were 
modeled into ITS CYLTRAN to compare its accuracy with the dosimeter readings. 
ITS CYLTRAN results are in agreement with the TLD measurements as discussed in 
the following chapter. 

CR-39 (LiF) neutron dosimeters were taped to the magnet sample-holder to measure 
the neutron dose prior to the scheduled irradiation periods. This vielded an average 
measured dose of 0.5 +0.3 megarads per coulomb in the sample plane [Ref. 9: p. 3]. A 


neutron dosimeter was also used to measure the neutron dose at the reference magnet 


I] 


position. According to its reading of 0.15 megarads per coulomb, the reference magnet, 
which had been mounted for the duration of the entire experiment on the same plane 
as the sample wheel holder but shielded from the bremsstrahlung radiation cone by 5 cm 
thick lead bricks, received a total neutron dose of 3.7 megarads, significantly less than 


gamma radiation exposure. 


B. EQUIPMENT 

The experiment conducted at LLNL exposed a total of 21 permanent magnets. 
Samples from four different manufacturers were irradiated. Table 2 lists the types of 
rare earth magnets used. Each sample had a cylindrical shape measuring 8 mm in di- 
ameter and 12 mm in length. The samples were mounted on the spokes of a water- 
cooled aluminum wheel pictured in Figure 4 on page 14. This sample wheel holder was 
part of a fully automated “in-situ magnetometer” apparatus shown in Figure 5 on page 


15 which was designed to remotely measure the remanence of each sample [Ref. 3: p. 3]. 


Table 2. REPM SAMPLES 


ALLOY MANUFACTURER B, (kilogauss) 


INGeie na Crucible Materials 4.85 
Wd as Hitachi Magnetics 4.78 
Sm, Coe Shin-etsu 4.35 


Preleainy Advanced Material Corp 4.58 





The magnetometer was designed to hold eight samples, but during the irradiation runs, 
only six holders were used. The reason for this was to achieve a significant dosage dif- 
ferential between the inner and outer three magnets. The two spokes with single sample 
holders did not provide adequate distance between the other magnets and therefore were 
not used. Each sample was water-cooled to eitiminate thermal effects. The temperature 
of each sample was monitored with resistive temperature detectors (RTD’s). Temper- 
ature and remanence measurements were interrogated periodically with a computer 


interfaced with the in-situ magnetometer apparatus. This was made possible by two 


pneumatic cylinders, controlled by relays and a stepping motor. See Figure 5 for the 
location of these control devices. A servo motor was designed to drive the coil platter. 
The coil platter (a lucite wheel resembling a film reel) moved between three discrete radii 
on the sample wheel platter by the pneumatic cylinders or positioners. The positioners 
were controlled by the computer. The sample platter was designed to rotate + 90° with 
the stepper motor. The individual samples were brought into a position where they 
could be interrogated by a pair of search coils (rotating coils) on the coil platter. The 
platter rotated freely (slowly) when the magnets were being irradiated. When magnets 
were being interrogated, the beam was secured off and the wheel rotation rate was in- 
creased to 7/S5 RPMs. The induced signals were brought out through low-noise, mer- 
curv slip rings to a gated integrator, digitized and recorded on computer disk. [Ref. 10: 
p. 2] 

A calibration to determine the accuracy of the in-situ apparatus measurement was 
performed with a set of samples that had been magnetized to different levels and meas- 
ured with the integrating fluxmeter. The calibration curve is shown in Figure 6 on page 
16. The one to two percent errors are due to small differences in the position of the same 
magnet with respect to the search coil pair. For each magnet sample, the repeatability 
between successive measuremenis was also on the order of one to two percent. The re- 
maining error could be due to dynamic misalignments between the vibrating sampling 
Wheel and the magnet samples. 

The charge deposited on the production target was measured directly with an Ortec 
Model 439 digital current integrator interfaced to an Ortec scalar to provide charge in 
coulombs. A conversion factor was applied based on the ratio of the average current 
measured at the production target to the average current measured at the beam dump 
without the production target, because the production target was less than range thick. 
meer. 9: p. 2] 

The open circuit remanence of each sample before and after irradiation was meas- 


ured using an integrating fluxmeter (Walker Scientific Model MF-3D). [Ref. 9: p. I] 


C. IRRADIATION RUNS 
1. Cobalt 60 Source 
Two samples of Sm-Co alloy were placed in a Co® tank to accumulate gamma 
radiation. The samples were Shin-etsu R26H ( Sm,Co,, ) and Vacomax 225 (SmCo,). 


Both samples were exposed to a total dose of 1.4 gigarads. Initial and post irradiation 
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Figure 6. Calibration of the In-Situ Magnetometer 


remanence readings were measured bv a Hall probe gaussmeter. Both samples experi- 
enced no experimentally resolvable remanence loss. 
One sample of Nd-Fe-B alloy (Crumax 355) was exposed to a total dose of 69 
megarads by the Co® source. No remanence loss was detected in this sample as well. 
The uncertainties in measuring the remanence of cach sample are due to the 
Hall probe gaussmeter sensor being comparable in size to the diameter of the magnets. 
Thercfore, the measurement is verv dependent upon the exact positioning of the probe 


relative to the magnet. Also, the measurement was observed to vary with temperature. 


2. Direct Electron Beam 

Two magnet samples were irradiated with a direct 82 MeV electron beam. The 
samples tested were Vacomax 225 (SmCo,) and Crumax 355 (Nd,Fe,,B). The beam pa- 
rameters used were set at an average current of 600 microamps using pulses of 2.8 
microsecond duration. A total of 25 pulses were deposited on the magnet samples for 
an estimated total dose of 36 kilorads. The Vacomax 225 measured no remanence loss 
whereas the Crumax 355 sample measured a loss of 1.5 percent in remanence field, which 
was a Statistically significant drop. This seems to suggest that a direct high energy 
electron beam causes more remanence loss than Co® gamma fields for the same total 
dose. 

The purpose of these measurements was to determine how much beam current 
would be required to induce a measurable loss 1n remanence for each magnet. This in 
turn would establish an estimated interrogation time interval for the irradiation runs 
utilizing the in-situ magnetometer apparatus. 

3. Bremsstrahlung Radiation 
An §5 MeV electron energy beam was used for all irradiation runs. Beam and 


other pertinent parameters are listed in Table 3 for the two scheduled periods. 


Table 3. LINAC OPERATING PARAMETERS 


26 JAN 88 16-19 FEB 88 


Distance from target (cm): 49.2 78.7 


Measured current at production target (uamps): 28.0 70.0 


Measured current at beam stop (uamps): 60.0 15.0 
Production target/beam stop current ratio: 2.14 Migs, 
Pulse repetition rate (Hz): 120.0 


Estimated beam spot diameter (cm): : 1.0 





Six samples of Hicorex 94EB (Nd,Fe,,B) were irradiated during run 1. (Refer 


to Table 4 on page 18 for a description of each irradiation run). An estimated total dose 
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of 0.4 gigarads was deposited on the magnet sample that was positioned on the beam 
axis. Post irradiation measurements on the integrating fluxmeter recorded a three to five 


percent drop in remanence. [Ref. 9: p. 4] 


Table 4. SUMMARY OF THE IRRADIATION RUNS 


Time Exposed Charge Number 
Run Nr Date Samples (hours) (coul) Irradiated 


26 Jan 88 Hicorex 94EB 3.0 2.0 


167 Feb So ee rumdneses 


18 FEB 88 Shin-etsu R26H 
AMC 


19 Feb 88 Hicorex 94EB 
AMC 


18/19 Feb 88  Shin-etsu R26H 
AMC 





Six samples of Crumax 355 (Nd,Fe,,B) were irradiated during run 2. When 
measurements on these samples were recorded on the second day, the interrogation 
readings measured zero remanence due to mechanical failure of the in-situ 
magnetometer. On the following morning, we discovered the lucite wheel containing the 
search coils cracked due to radiation damage. The damaged wheel caused erroneous 
zero remanence measurements. After measuring the post irradiation remanence of these 
samples on the integrating fluxmeter, onlv a two percent loss was measured. A calcu- 
lated dose of 0.92 gigarads was deposited [Ref. 9: p. 4]. 

During run 3, the six Crumax 355 samples were replaced by four samples of 


Shin-etsu R261] (Sm,Co,,) and two samples of AMC (Pr,,Fe-,B,). These samples were 


irradiated for a total of seven hours. Active cooling of the samples was no longer fea- 
sible since the water circulation svstem also suffered radiation damage. 

Two Shin-etsu R26H samples and one of the AMC samples were replaced with 
two samples of Hicorex 94EB and a different sample of AMC during runs 4 and 5. The 
remaining two samples of Shin-etsu R26H and the other AMC sample remained in their 
respective positions to accumulate more dose. The two Shin-etsu R26H and the AMC 
samples that were irradiated for two days (15 hours) received a total dose of 1.5 and 1.8 
gigarads respectively. Post irradiation measurements on the integrating fluxmeter 
measured no loss in remanence for the Shin-etsu R26H sample. The AMC samples re- 
ceived the highest dose total and recorded the biggest drop, a 75 percent loss in 


memenence. [Ref. 9: p. 5) 


D. TEMPERATURE EFFECTS 

For the first three scheduled days of irradiation runs, the temperature of each sample 
was monitored. Temperature recorded for each sample was maintained at around room 
temperature. Therefore. remanence loss for the Hicorex 94EB and Crumax 355 samples 
Were not due to temperature. For the irradiation runs conducted on February 18 and 
19, the water cooling system was not available and temperature was not monitored. The 
reference magnet Was not water-cooled during the entire irradiation period. When the 
in-situ magnetometer apparatus was operational, an average temperature of 26 °C was 
measured after four hours of irradiation. Based on this information, an estimated tem- 
perature for the samples irradiated on February 18 and 19, after the apparatus was 
damaged. reached a temperature slightly higher than 26 °C, since total irradiation time 
Was greater, cight hours, and the magnet samples were exposed directly to the 
bremsstrahlung radiation. The temperature still was maintained well below the stabiliz- 
ing temperature of each magnet sample. Therefore, temperature had minimal effect 


throughout the entire experiment. 


Ey RESULTS 

Table 5 on page 20 summarizes the irradiation results. Crumax 355 (Nd,Fe,,B) was 
irradiated using a direct electron beam. A loss in remanence was measured which could 
have been due to the positioning error of the hall probe. The values for the normalized 
remanence at maximum dosage for the samples exposed to bremsstrahlung radiation 
were computed using the following linear fit equations used to plot the curves in 


Fie@re 7 on page 21: 
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where a = -0.14 for Hicorex 94EB)=] 25 x 1072 for Crumax 355, 2.25 x@l@> for Shin-eetse 
R26H and -0.33 for AMC. 


Table 5. RESULTS OF IONIZING RADIATION ON REPM SAMPLES: Val- 
ues for bremsstrahlung radiation corresponds to the fit values from 
Figure 7 on page 21. 


Normalized 
Temp Maximum Dose  Remanence 
Sample Alloy (°C) (10° rads) at Max Dose 
Cobalt 60 Source: 
VACO e725} smCo, 24 1.400 1.000 
SHANE Sia cr Sm Co. 24 1.400 1.000 
CROMAX 395 PGs cb 24 069 1.000 
Direct Electron Beam: 
VY ACONTA™ 225 5m,Co,- nm 36.0krads 1.000 
CRUST AY 355 ial Srey 8: nm 36.0krads 985 
Bremsstrahlung Radiation: 
HICOREN 94EB dics 24 -396 945 
CRE Nee 355 NCH alee: 23 ols OU 
SLING PS Ge Oil 5m,C0,; 28 Pos 1.000 
AMC Pijigessny 28 1.782 412 


nin - not measured 


Post irradiation remanence was measured on an integrating fluxmeter and normal- 


ized. ITS,CLYTRAN simulation codes were used to plot the normalized remanence 
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Figure 7. ITS/CYLTRAN Plot of Normalized Remanence vs Dose 


versus dose in gigarads for each magnet alloy. Energy deposition in gigarads per 
coulomb was computed using CYLTRAN and converted to gigarads based on the total 
charge deposited on each magnet sample recorded during the irradiation runs. A linear 
fit as previously described was used to plot the curves for each magnet sample. 

The same procedure was used to plot the results obtained from the irradiation runs. 
The difference between the two plots was that the total dose computed was based on the 
flash x-ray calculation described in the Preliminary Calibration Measurements section 


of this chapter. Figure 8 on page 23 represents the plot which lays out the general trend 


of remanence loss for each magnet sample irradiated and dose calculated using the flash 
x-ray scaling. A linear fit was also used to plot the curves for each magnet sample with 
the exception of the AMC sample. Refer to Appendix C for a description on how 
ITS'CYLTRAN codes are used and executed. 

Figure 8, which is the plot produced from the flash x-ray calculation is very similar 
to Figure 7, the plot produced from the ITS/CYLTRAN simulation. The only signif- 
icant difference that exist 1s the nonlinear curve for the AMC sample shown in 
Figure 8. The disparity could have been attributed to the uncertainties in the flash x-ray 
calculation used to compute the total dose. 

In general, ITS’CYLTRAN computed higher on-axis dose for each magnet alloy 
compared with the data obtained from the flash x-ray calculation. I[TS/CYLTRAN, 
however, computed a slightly less dose for those magnets positioned on the outer radius 
of the magnet sample holder. (See Figure 9 on page 24 for the magnet positions on the 
sample holder). Table 6 on page 26 shows the difference in dosage values computed by 
ITS/CYLTRAN and the flash x-ray scaling used at LLNL for the different magnet 
samples positioned on the sample holder. 

For the samples that suffered noticeable demagnetization, a change in the coercivity 
field might occur. This can be examined by plotting a demagnetization curve in a 
hysteresis graph. Future studies in terms of the Barkhausen effect [Ref. 11: p. 46], which 
is an experimental characterization of the dvnamics of domain wall movement, precludes 
demagnetizing the irradiated magnets to obtain the hysteresis loop required to determine 
the coercive force of the sample. Studies of the Barkhausen effect on the irradiated 
samples are planned for June and July, 1985 to investigate the causes to why radiation 
exposure reduces the remanence i pemmancnt maeneisajicie pall 

A Crumax 355 (Nd-Fe-B allov) was used and positioned along the same plane with 
the sample-holder. This reference magnet was shielded from the bremsstrahlung cone, 
but still received one-third to one-half the neutron dose of the other samples since neu- 
tron radiation is not as directional. Exposure to neutron dose was the only means of 
remanence loss in which the reference magnet could have encountered, since the refer- 
ence magnet was shielded from bremsstrahlung radiation. Since post irradiation meas- 
urements recorded no remanence loss, we can conclude that neutron fluence had no 
significant effect on remanence loss on the magnet samples irradiated based on the 


dosimeter measurements taken prior to irradiation. [Refi 9: p. 6] 
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Figure 8. Flash X-ray Scaling Plot of Normalized Remanence vs Dose 


F. CONCLUSIONS 

In conclusion, the experimental results show that the Shin-etsu R26H (Sm,Co,-,) 
magnets are the most radiation resistant of the materials studied. Sm,Co,. provides the 
greatest hardness against bremsstrahlung irradiation. ITS'CYLTRAN analvsis in the 
next secuon will show that gamma radiauon is the dominant field produced from the 
production target. 

Comparing the results of the two Nd,Fe,,B alloys, Crumax 355 and Hicorex 94EB, 


Hicorex 94EB proved to be less resistant. The magnets were manufactured by different 





Figure 9. Niagnet Positions on the Sample Wheel Holder: where the radius of the 
inner Magnets 1s 2.54 cm and the outer magnet is 6.55 cm. 


companies, Which confirms that manufacturing techniques of similar allovs have impact 
on radiation hardness. Barkhausen effect studies may provide added insight to these 
differences. 

Based on the neutron dosimeter readings placed at each sample site, including the 
reference magnet, remanence loss due to neutron exposure was insignificant. 

The temperature effect on permanent magnets was held minimal to none, since the 


Curie temperature for each magnet alloy is far above its operating temperature during 


radiation exposure. The Curie temperatures for iron and cobalt are 770 °C and 1131 
“Grrespectively’. 

During the planning phase, prior to the irradiation runs, numerous test and equip- 
ment (1.e. in-situ magnetometer interfaced with computer system) calibrations were made 
to ensure accurate data would be recorded throughout the experiment. We continued 
the experiment in spite of equipment failure and kev results and conclusions were ob- 
tained during these final davs of scheduled irradiation runs. Our results indicate that to 
acquire the gross radiation effect measurements, the sophisticated apparatus used was 
not needed. However, it would be most useful in studies where prototvpe magnets of 
different compositions are simultaneously exposed. This would allow us to mix different 
magnet samples on the same wheel and expose them to the same radiation environment 


and trace their demagnetization curves out to significant remanence losses. 
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IV. DOSIMETRY ANALYSIS USING ITS/CYLTRAN 


ITS CYLTRAN Monte Carlo simulation codes were used to describe the pro- 
duction and transport of electrons and photons through materials of different materials 
of cylindrical geometry. These codes were utilized to determine electron and photon 
energy spectra, angular spectra and resulting energy deposition in materials [Ref. 12: p. 
6]. Appendix C describes the use and execution of ITS/CYLTRAWN codes. 

The experimental configuration at LLNL’s 100 MeV Linac was modeled with 
ITS, CYLTRAN codes to calculate the dose deposited on the three different magnet al- 
loys. A comparison was made between the dose calculated using CYLTRAN and the 


measured results obtained with dosimeters. The result compare favorably. 


A. ENERGY DEPOSITION 

CYLTRAN was used to calculate energy deposition on a circular disc one centime- 
ter thick and 10 cm in radius made of a given magnetic material. Sm,Co,,, Nd,Fe,,B and 
Pr,,Fe-,B, were analvzed during different runs. The radius corresponded to the radius 
of the sample wheel holder used in the experiment. Ten subzones were used to calculate 
the dose at one centimeter intervals from the center of the disc. This number of sub- 
zones allowed the dose to be determined at the corresponding magnet position on the 
sample holder. (Refer back to Figure 9 on page 24 for the magnet positions on the 
sample holder). 

Angular divergence and beam energy smearing were common operating character- 
istics of the beam at LLNL’s Linac and were modeled into the simulation. Detailed 
descriptions of angular divergence and beam energy smearing will be given later in this 
section. An experimentally determined beam radius 0.5 cm was used in all 
ITS’'CYLTRAN simulation runs. Table 7 on page 34 through Table 9 on page 38 give 
the computed doses of primarv electrons, photons and the total energy deposited used 
to plot the radiation dosage patterns. Figure 10 on page 35 through Figure 12 on page 
39 provide the radiation dosage patterns for the magnet samples. (All tables and figures 
referenced are at the end of this chapter). The total energy deposited in each magnet 
includes contributions from primary and secondary electrons as well as gamma radi- 
ation. The number next to each dose value represents the calculated uncertainty. It 


represents a measured percentage of the calculated value. ITS calculates the uncertainty 


Zi 


based on stochastic processing as each incident electron traverses a segment of the me- 
dium. This process is initiated based on the number of histories and batches simulated. 

As shown in Figure 10 on page 35 through Figure !2 on page 39, the shape of the 
dosage patterns are relatively the same. The slight difference in values can be attributed 
to the difference in density of each alloy. The measurcd densities forssnne@o- 
KNd,Fe,,B and Pr,,Fe:,B, are 8.395 7:60)-and 10°20 srcmie respectively. Sip@a mae 
Pr,.Fe..B,, which has a higher density, received a slightly higher dose compared with 
Nd,Fe,,B. Notice also that for each alloy, bremsstrahlung radiation is the dominant 
field. Since tungsten has a high atomic number, the bremsstrahlung radiation is corre- 
spondingly high. The fraction of energy loss of primary electrons that are given off as 
photons are more dominate as electron energy increases. This energy loss by radiation, 
known as bremsstrahlung, 1s proportional to Z? where Z is the atomic number material 
and increases linearly with the electron energy. This accounts for the relatively low en- 


ergy depositcd by primary electrons as shown on each graph. 


B. DOSAGE DEPOSITION ON VARIOUS MATERIALS 

Energy deposition on various materials was modeled using ITS/CYLTRAWN to de- 
termine how the absorbed dose at the beam axis varies in different materials. The sim- 
ulated radiation dosage was produced by an 85 MeV electron beam of radius 0.5 cm on 
a production target similar to the one used at LLNL’s Linac. Table 10 on page 40 list 
the types of materials and alloys used in these simulation runs and compares them with 
the relative dosc to silicon. The comparison with silicon was made because silicon is a 
common elcment used in dosimeters to measure radiation effects. Silicon is also the 
material used in semi-conductors, which is one of the materials of primary interest with 
respect to radiation hardness. Using the same operating parameters for each material, 
the absorbed dose deposited on the common material used in dosimetry, Ca, Si and LiF 


are a factor of three less than those elemcnts and alloys used to manufacture magnets. 


C. PRELIMINARY MODELING 

Prior to establishing the final model just previously discussed, a series of preliminary 
runs were conducted to compare the differences between various input parameters. The 
first comparison made was between the electron energy spectra of tungsten and the 
conmiposite material “heavimet” (4% Cu, 7% Ni, 89% W). Since the production target 
used at LLNL’s Linac was heavimet, the energy and angular spectra between the 
heavimet and pure tungsten were generated to determine any differences. Figure 13 on 


page 4] compares the cnergy spectra gcnerated through a tungsten disc and heavimet 
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using a point source energy beam at 85MeV. More electrons were transmitted in the 
lower energy spectrum of the tungsten foil compared to heavimet, leading us to the 
conclusion that tungsten produces more low energy secondaries. Looking at the oppo- 
site end of the spectrum, heavimet allowed a greater number of high energy electrons 
(primaries). These results were expected based on the density of each target. Heavimet 
has a density of 17.1 gr/cm*’ compared to tungsten at 19.3 gr/cm. The radiation length 
of a material depends upon its density (for constant atomic number). Since both mate- 
rials had the same thickness, and the density of heavimet was slightly less than tungsten, 
more primary electrons would be transmitted through heavimet. 

Another comparison run was made to determine what effects beam radius has on 
the energy spectrum of the heavimet disc. A point source and a monoenergetic electron 
beam of radius 0.5 cm were compared. The energy and angular distributions (see 
Figure 14 on page 42 and Figure 17 on page 45), exhibited no significant differences, 
Which implies that beam spot size has minimal effect on the energy and angular spectra 
for the eiectrons and photons. This is obviously true onlv when the beam radius is less 
than the radius of the heavimet disc. Later in this section, an analysis of beam spot size 
will be conducted on the actual production target. A difference in primary electrons and 


gamma radiation dosage patterns will be observed. 


D. ANGULAR DIVERGENCE 

Angular divergence considers the angular distribution of the source particles with 
respect to the beam axis. The divergence angle was determined by measuring the effec- 
tive aperture of the quadrupole magnets and the diameter of the production target, as 
Well as the distance between them. An angle of 0.18° was calculated based on these di- 
mensions. Angular divergence was modeled by CYLTRAN. Using a 0.5 cm beam ra- 
dius and an 85 MeV electron energy beam for the electron source, Figure 18 on page 
46 through Figure 21 on page 49 compares the energy and angular energy spectra for 
transmitted electrons and photons through heavimet for a beam with and without an- 
gular divergence. No significant difference was noted which was expected since the d1- 


vergence angle was significantly small. 


E. ENERGY SMEARING 

Energy smearing is defined as the energy distribution describing the electron beam. 
Energy smearing at LLNL’s Linac resembles a Gaussian distribution with a five percent 
energy drop occurs at the full width half maximum (FWHM) point on the Gaussian 


distribution curve. A five percent energy smearing was modeled into CYLTRAN to 


yay 


observe the effects it would have on the energy spectrum for transmitted electrons. 
Continuing to use a 0.5 cm beam radius and an 85 MeV electron energy beam for the 
electron source, Figure 22 on page 50 through Figure 25 on page 53 compares the en- 
ergy and angular spectrums for transmitted electrons and photons through heavimet for 
a beam with and without energy smearing. The onlv significant difference noted 1s the 
gradual decrease of the number of electrons/MeV at the higher energy bins for the en- 
ergy distribution modeled with energy smearing. This was expected since the energy at 
the outer diameter of the electron beam decreases to 79 MeV. Since a noticeable dif- 


ference occurred, energy smearing was modeled into the remaining CYLTRAWN runs. 


F. HEAVIMET DISC VS THE PRODUCTION TARGET 

A full scale drawing of the longitudinal cross section of the production target used 
during the irradiation runs 1s shown in Figure 2 on page 10. The 2.5 mm thick heavimet 
disc is cooled by an 8 mm thick copper laver support structure. A comparison run was 
simulated using the heavimet disc with and without the copper layer to determine how 
the energy and angular spectrums differed for electrons and photons. Angular diver- 
gence and energy smearing was included into the modeling. Figure 26 on page 54 
compares the energv spectrum for the heavimet disc and the production target with the 
copper laver. As shown, a significant difference exist at the higher energy bins. A de- 
crease in the number of electrons’ MeV occurred when the copper laver was added. This 
suggest that the combined thickness of the heavimet disc and the copper laver made it 
more difficult for the high energy primary electrons to be transmitted. The majority of 
the transmitted electrons were the low energy secondaries. By adding the 8 mm copper 
laver to the heavimet disc, the number of radiation length of the material was increased 
by over 40 percent. The electron energy loss due to radiation was calculated for both 
cases. The calculated electron energy loss for a pure heavimet material was 6.1 MeV 
compared to 16.4 MeV for the production target [Ref. 13: p. 68]. Refering to Figure 26 
on page 54, notice the curve representing the production target loses all its energy at 
around 69 MeV, and the curve for heavimet loses energy significantly at about 74 MeV. 
This corresponds almost exactly with our calculations. This gives us confidence in the 
results generated by ITS’CYLTRAN codes. Figure 27 on page 55 through Figure 29 
on page 57 are the remaining comparisons between energy and angular spectrums for 
electrons and photons. Figure 27 on page 55 and Figure 29 on page 57 compare the 
angular distributions. Note what happens when the copper support structure is added 


to the heavimet disc. The half-width angle increases for both electron and photon 


angular distributions. This 1s because the radiation length for copper is almost twice as 
much as tungsten, 12.86 gr/cm* versus 6.76 gr/cm? respectively, Which accounts for the 
greater angular spread of transmitted photons and electrons. Radiation length 1s 
inversely proportional to the deflection angle of transmitted photons and electrons which 
accounts for the increase in half-width angle. 

The half-width angle corresponds to the angle for which the 1/e (.368) of the maxi- 
mum number of electrons or photons per stere-radian per MeV is at. This value can be 
utilized to position any equipment, apparatus or any material from the radiation source 
in order to receive a desired dose. Preliminary modeling using ITS/;CYLTRAN can 


greatly aid an experimental setup to achieve optimum results. 


G. HALO EFFECT 

The halo effect is when the electron beam cannot be focused to a size smaller than 
the radius of the heavimet disc. As a result, the beam initially hits the copper support 
structure before the heavimet foil. producing a halo of electrons. Since the radiation 
length of copper is greater than tungsten, more electrons will transmit through the 
copper material. A higher exposed dose of electrons will occur at the outer radius ofa 
sample material. This was evident during one of the preliminary experiments prior ‘to 
the irradiation runs. Glass plates were known to darken when exposed to high levels of 
electrons. The electron beam during the preliminary experiments were estimated to be 
greater than the radius of the heavimet foil, and as a result, darken rings were deposited 
on the glass plates. An electron beam radius of 1.00 cm was modeled into CYLTRAN 
to verify the occurence of a high level of electron dose. The radius of the heavimet disc 
is 0.943 cm, which is less than the radius of the simulated electron beam. A simulated 
run was executed to measure the dose on a disc platter made of LiF (a common material 
used in TLD chips to measure gamma dose). The radius of the LiF platter was modeled 
at 10 cm to correspond with the radius of the sample wheel holder. Figure 30 on page 
58 shows the dose pattern for the primary electron, gamma radiation and total energy 
deposited in gigarads per coulomb. As expected, the primary electrons dominate at the 
outer radius. Figure 31 on page 59 shows the radiation dosage pattern for an energy 
beam with radius 0.75 cm. Note that the primary electrons dominate only at the outer 


perimeter of the LiF platter. 


H. EXPERIMENTAL VERIFICATION 
The only experimental result which we can use to compare ITS/CYLTRAN with 


was the dose measurement conducted prior to one of the irradiation runs. An array of 


S| 


LiF TLD chips were positioned 167.8 cm from the production target. The dosimeters 
measured 22 +5 megarads per coulomb along the beam axis. This same arrangement 
was modeled using the ITS/CYLTRAN simulation codes. A 20.4 megarads per coulomb 
On axis dose was calculated which concurs with the TLD readings. Figure 32 on page 
60 shows the radiation dosage pattern for LiF as it pertains to this configuration. The 


corresponding dosage pattern values are listed in Table 11 on page 61. 


I. DISCUSSION OF ERRORS 

The calibration curve shown back in Figure 6 on page 16 was used to measure the 
accuracy of the in-situ magnetometer apparatus. The curve shown was plotted using a 
linear fit. The one to two percent error was probably due to small differences in the 
positioning of the sample magnets with respect to the search coils. A part of the error 
reflects the 0.4 percent error in maintaining a precise rotation rate for the sampling 
wheel. The remaining error could be due to dynamic misalignments between the vi- 
brating sampling wheel and the samples. [Ref. 9: p. 2] 

The uncertainty in the dosimetry calculations computed by ITS/CYLTRAN runs 
was approximately 10 percent of the calculated dose. The uncertainty due to the mod- 
eling of the dose pattern profile is difficult to quantify. The spread in the data shown 
back in Figure 7 on page 21 may reflect the overall uncertainty, however, we cannot rule 
out the sample to sample difference in magnets in general, regardless of the manufac- 


turing techniques and the manufacturer. 


J. CONCLUSIONS 
ITS;CYLTRAN has proven to be a useful tool in analyzing the energy and angular 
spectra for varying electron beam parameters. The following conclusions can be made 


from the series of angular and energy spectrums for the electron and photon plots 
produced by ITS/CYLTRAN: 


Beam spot size: No significant difference was noted between a point source energy 
beam and a monoenergetic beam of radius 0.50 cm. (See Figure 14 through 


Figure 17). This had no effect on the on-axis absorbed dose on a material. 


Angular divergence: There were no differences noted in the energy and angular spectra 
of electrons and photons, therefore there were no effect on theon-axis ab- 


sorbed dose on a material. (See Figure 18 through Figure 21). 


Energy smearing: A noticeable difference in the electron energy spectrum was noted 
when energy smearing was modeled into the simulation. A gradual decrease 
in the number of electrons per MeV occurred at the higher energy bins for 
the energy distribution. (See Figure 22 through Figure 25). Energy 
smearing also had an effect on the dosage pattern of a material. An 18 
percent decrease in the absorbed dose was computed when energy smearing 


was included. 


Heavimet Disc vs Production Target: The production target is composed of the com- 
posite material heavimet and copper. By adding the copper layer to the 
simulation, a significant decrease in energy occurs at the higher energy bins 
and increases at the lower energy bins. The half-width angle for the electron 
angular spectrum also increased when copper was added. (See Figure 26 


through Figure 29). 


We have verified the validity of ITS; CYLTRAN simulation codes with great accu- 
racy by dosimetry measurements. Therefore the dose measurements calculated by 
ITS'CYLTRAN are considered to be more accurate than the data computed using the 


flash x-ray calculation. 


oe 


Table 7. RADIATION DOSAGE VALUES AT SM-CO ALLOY: These values 
are used to plot the radiation dosage pattern for Figure 10 on page 35. 
Numbers next to each dosage value represents the measure of uncertainty 
in percent. 


DISTANCE PRIMARIES GAMMAS TOTAL 
(CM) (GRads/Coul) (GRads/Coul) (GRads/Coul) 


0.5000 0.0514 41 0.2613 10 0.3125 11 
1.5000 0.0149 26 0.1861 6 0.1988 7 
2.5000 0.0186 19 0.1309 0.1480 
3.5000 0.0244 17 0.1282 
4.5000 0.0229 14 0.1072 

O00 0.0194 14 0.0869 

6.5000 0.0174 10 0.0742 
7.5000 0.0204 13 0.0626 
8.5000 0.0230 8 0.0600 - 
9.5000 0.0229 11 0.0476 
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Figure 10. Radiation Dosage Pattern at Sm-Co Alloy: The dosage pattern was 
calculated using an 85 MeV electron beam with a radius of 0.75 cm. 
The absorbing material was posiuoned 78.7 cm from the production 
target based on 10,000 histories. 
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Tabie 8. RADIATION DOSAGE VALUES AT ND-FE-B ALLOY: Same de- 
scription from Table 7 applies accept values are used to plot the radiation 


dosage pattern for Figure 11 on page 37. 


DISTANCE PRIMARIES GAMMAS TOTAL 
(CM) (GRads/Coul) (GRads/Coul) (GRads/Coul) 
0.5000 0.0200 46 0.2311 1 0.2502 12 
1.5000 050275025 0.1600 10 0.1880 9 
2.5000 0.0241 22 0.1368 9 0.1603 8 
3.5000 020239 19 0.1097 6 UA ope aes) 
4.5000 0.0199 19 0.0937 7 0.1124 8 
5.5000 0.0277 10 0.0821 7 0.1093 6 
6.5900 0.0207 18 0.0717 8 Onn 2 7 
7.5000 0.0206 10 0.0585 6 0.0783 4 
8.5000 0.0162 11 0.0322 6 0.0677 6 
9.5000 0.0202 7 0.0460 7 UrU05/ 5 
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Figure 11. Radiation Dosage Pattern at Nd-Fe-B Alloy: Same operating parame- 
ters as Figure 10 on page 35. 
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Table 9. RADIATION DOSAGE VALUES AT PR-FE-B ALLOY: Same de- 
scription from Table 7 on page 34 applies accept values are used to plot 
the radiation dosage pattern for Figure 12 on page 39. 


DISTANCE PRIMARIES | GAMMAS TOTAL 
(CM) (GRads/Coul) (GRads/Coul) | (GRads/Coul) 


0.5000 0.0313 27 0.2267 8 
1.5000 0.0123 36 0.1874 10 
2.5000 0.0139 24 + 
3.5000 0.0167 15 
4.5000 0.0250 17 
5.5000 0.0189 16 
6.5000 0.0199 9 
7.5000 0.0234 8 
8.5000 0.0203 I1 
9.5000 0.0181 11 
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Figure 12. Radiation Dosage Pattern at Pr-Fe-B Alloy: Same operating parame- 
ters as Figure 10 on page 35. 
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Table 10. DOSAGE ON VARIOUS MATERIALS 


Density On-axis Dose Dose Relative 
Materia! (gr/cm?) (x 10° rads/coul) to Silicon 


Ca 55 0.104 0.97 
Si 0.107 1.00 
0.084 0.79 
0.352 
0.347 
0.300 
Sm,Co;; 3 0.289 
Nd,Fe,,B 0.250 
Pr,.Fe-5B, 0.285 
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Figure 13. Tungsten vs Heavimet Electron Energy Spectrum: The spectrum was 
calculated using an 85 MeV electron beam with a radius of 0.50 cm. 
The absorbing material was positioned 78.7 cm from the production 
target based on 10,000 histories. 
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Figure 14. Point Source vs Mionoenergetic Beam Electron Energy 


Spectrum: Same operating parameters as Figure 13 on page 41 used 
tO compare a point source beam with a monoenergetic beam of radius 
0.50 cm. 
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Spectrum: Same operating parameters as Figure 14 on page 42. 
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Figure 16. Point Source vs Monoenergetic Beam Photon Energy Spectrum: Same 
opcrating parameters as Figure I4 on page 42. 
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Figure 17. Point Source vs Mionoenergetic Beam Photon Angular 


Spectrum: Same operating parameters as Figure 14 on page 42. 
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Figure 18. Angular Divergence Comparison for Electron Energy Spectrum: Same 
operating parameters as Figure 14 on page 42 accept a divergence an- 
gle of 0.18° was incorporated into the simulation. Comparison was 
made with a monoenergetic beam without angular divergence. 
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Figure 19. Aneular Divergence Comparison for Electron Angular Spectrum: Same 
operating parameters as Figure 18 on page 46. 
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Figure 20. Angular Divergence Comparison for Photon Energy Spectrum: Saime 
Operating parameters as Figure 18 on page 46. 
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Figure 21. Angular Divergence Comparison for Photon Angular Spectrum: Same 
operating parameters as Figure 18 on page 46. 
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Figure 22. Energy Smearing Comparison for Electron Energy Spectrum: Same 
operating parameters as Figure 18 on page 46 accept a divergence an- 
gie of 0.18° and five percent energy smearing was incorporated into the 
simulation. Comparison was made with a monoenergetic beam with 
angular divergence. 
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Figure 23. Energy Smearing Comparison for Electron Angular Spectrum: Same 
operating parameters as Figure 22 on page 50. 
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Energy Smearing Comparison for Photon Energy Spectrum: 


Operating parameters as Figure 
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Energy Smearing Comparison for Photon Angular Spectrum: 
operating parameters as Figure 22 on page 50. 
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Figure 26. Heavimet Disc vs Production Target Electron Energy Spectrum: Same 
Operating parameters as Figure 22 on page 50. Compares spectra with 
Q.60 cm copper lining added to the production target. 
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Spectrum: Same operating parameters as Figure 26 on page 54. 
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Figure 28. Heavimet Dise vs Production Target Photon Energy Spectrum: Same 
operating parameters as Figure 26 on page 54. 
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Heavimet Disc vs Production Target Photon Angular Spectrum: Same 
Operating parameters as Figure 26 on page 54. 
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Figure 30. Radiation Dosage Pattern for LiF with a Halo Effect: The dosage 
pattern was calculated using an 85 MeV electron beam with a radius 
of 1.00 cm. The absorbing material was positioned 78.7 cm from the 
production target based on 10,000 histories. 
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Figure 31. Radiation Dosage Pattern for LiF without a Halo Effect: Same oper- 
ating parameters as Figure 30 on page 58 accept beam radius of 0.75 
cm was simulated. 
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Figure 32. Verification of the Radiation Dosage Pattern for LiF: The dosage 
pattern was calculated using an §5 MeV electron beam with a radius 
of 0.75 cm. The absorbing material was positioned 167.8 cm from the 
production target based on 10,000 histories. 
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Table 11. RADIATION DOSAGE VALUES FOR LIF TLD CHIPS 


DISTANCE PRIMARIES GAMMAS 


(GRads/Coul) 


TOTAL 
(GRads/Coul) 


(CM) (GRads/Coul) 


0.5000 
1.5000 
2.5000 
3.5000 
4.5000 

-5000 
6.5000 
7.5000 
8.5000 


9.5000 


0.0002 90 
0.0015 61 
0.0014 44 
0.0009 34 
0.0017 18 
0.0012 27 
0.0015 27 
0.0014 31 
0.0013 22 


0.0014 32 
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0.0203 18 
0.0166 8 
0.0198 7 
0.0134 6 
0.0121 
0.0116 
0.0110 
0.0090 
0.008 1 


0.0061 


0.0204 18 
0.0180 8 
0.0213 5 
0.0143 6 
0.0138 - 
0.0128 

0.0125 

0.0103 10 
0.0095 7 
0.0075 8 





V. CONCLUSIONS 


After radiating a total of 21 permanent magnet samples from four different man- 
ufacturers and exposing them to bremsstrahlung doses of up to 1.8 gigarads, the fol- 
lowing conclusions can be made based on the irradiation runs conducted at LLNL’s 100 
MeV Linac: 


1. Sm,Co,, was found to be significantly less sensitive to gamma radiation than 


Nd,Fe,,B, which agrees with previous neutron and proton irradiation experiments. 


2. Temperature effects were insignificant and did not contribute to the magetization 
loss of any sample. For the initial irradiation runs, the temperature of each sample 
was maintained near room teniperature by a water-cooling system. An out-of- 
beam reference magnet which was not water cooled during the irradiation runs, 
measured no loss in remanence. The monitored temperature at the reference mag- 
net measured slightly above room temperature of 28 °C. We can therefore con- 


clude that temperature did not reach a level high enough to cause remanence losses. 


o>) 


The samples of Nd,Fe,,B were provided by two different manufacturers, Hitachi 
Magnetics Corporation ({llicorex 94EB) and Crucible Magnetics Company 
(Crumax 355). The Hicorex 94EB samples lost three to five percent of ther omp 
ginal renianence compared to onlv a one to two percent drop for the Crumax 355 
after absorbing a maximum dose of 0.92 gigarads. Samples from different man- 

facturers exhibit large differences in the demagnetization as a function of dose. 
Different production processes, the presence of different alloy additives or the 
microcrystalline structure could be responsible for the difference in the rate of 


remanence loss. 


4. The sample alloy of Pr,,Fe.,B, proved to be extremely sensitive to gamma radiation 
losing 75% of it origina! remanence after absorbing a maximum dose of 1.78 


gigarads, 


Future studies should further explore these trends. Development of new magnetic 
materials have provided a more promising outlook in magnet technology. The man- 


ufacturers of permanent magnets should continue to provide new material samples to 


radiation facilities to continue this field of research and determine the optimum perma- 
nent magnet for use in high radiation environments. As the data base expands, 
phenomenological models, which address microstructural factors may evolve to levels 
that provide guidelines for the manufacture of radiation resistant alloys. A comprehen- 
sive program to assess phenomenology of many types of magnets would benefit from the 
use of an automated remote “in-situ” magnetometer. 

The ITS analysis of the experimental data obtained from the irradiation runs proved 
to be extremely valuable. The comparisons between measured results and ITS simu- 
lation runs showed high levels of agreements. This proved that ITS can be utilized more 
extensively to analvze similar models and achieve useful results. With the operating cost 
of irradiation time steadily increasing, utilization of the ITS/CYLTRAN Monte Carlo 
simulation codes to analyze energy deposition, energy and angular distributions can 
prove extremely economical. More sophisticated simulation codes, which utilizes the 
same Monte Carlo scheme to calculate electron, photon and neutron spectra are in use 
today. By taking advantage of these resources, more in-depth analysis of radiation ef- 
fects can be conducted at significantly lower cost. 

Permanent magnet technology and their applications have grown to become a S$] 
billion industry. Though improvements in the quality and strength of magnetic materials 
continue, Sm,Co,- still remains the best suited magnet compound for high radiation en- 
vironment. 

Based on our recent study and the previous studies conducted at LANL, CERN and 
TRIUMPF, these radiation damage experiments lead us to the empirical observation 
that the order of radiation hardness is Sm,Co,,, SmCo, and Nd,Fe,,B, regardless of the 


source of radiation, 1.e. gammas, electrons, protons and neutrons. 


63 


APPENDIX A. MAGNETIC PROPERTIES OF MATERIALS 


The following concepts and terminologies of magnetic properties are provided to aid 
the reader in obtaining an introductory understanding of magnetic materials. The fol- 
lowing discussion has been adapted from Donald R. Askeland’s book titled, The Science 
and Engineering of Materials. Askeland’s description are concisely written for anyone 


who has a minimal background in the study of magnetic properties of materials. 


A. ATOMIC STRUCTURE OF MATERIALS 

Atoms comprised of a positively charged nucleus are surrounded by electrons which 
move about the nucleus in definite states of total angular momentum. Studies have 
shown that magnetiSni occuring in transition metals is associated almost entirely with 
the states of electrons in the third shell of the atomic structure. The ordering of spins 
of the electrons produces magnetic fields as if charges were in motion. 

In most materials there is an identical number of electrons with opposite orientation. 
Since the polarity of the magnetic field is determined bv the direction of the electron 
spin, equal but opposite spins cancel the inagnetic field. Thus no magnetic field exist. 

Ferromagnetic materials have high levels of magnetization due to unpaired 
electrons. Associated with each electron is a magnetic moment. This moment, called 


the Bohr magneton 1s equa! to [Ref. 14: p. 607] 


gh 
dim 





~24 2 
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where q 1s the charge on the electron, h is Planck’s constant and m, is the mass of an 
electron. In an atom, the magnetic moment 1s due to the orbital angular momentum 
of the electron as well as the spin of the electron (Figure 33 on page 65). [Ref. 14: p. 
608] 

Atoms consist of a specific electronic structure characterized by definite quantum 
state and quantum numbers. Each quantum state can contain two electrons, each hav- 
ing an opposite spin. The magnetic moments of each electron pair 1 a quantum state 
are opposed and, consequently, whenever a state is completely full, there is no net 
magnetic moment. In contrast, any atom of an element with an odd atomic number can 
have a net magnetic moment from the unpaired electron. However, this is generally not 


true. In most of these clements, the unpaired electron is a valence electron. Because the 
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Figure 33. Origin of Magnetic Dipoles: (a) The spin of the electron produces a 
magnetic field with a direction dependent on the quantum number m,. 
(b) electrons orbiting about the nucleus create a magnetic field about 
the atom. 


Valence electrons from each atom interact in forming the bulk material, the magnetic 
moments, statistically cancel and no net magnetic moment is associated with the mate- 
rial. Certain elements, such as the transition metals, have an inner energy level that 1s 
not completely filled. Figure 34 on page 66 shows the electron spins in the 3d energy 
level in transition metals for the elements scandium through copper. The arrows indicate 
the direction of spin. [Ref. 14: p. 609] 

Tne electrons in the 3d level of the remaining transition elements do not enter the 
shell in pairs. Instead. as in magnanese, the first five electrons have the same spin. Only 
after half of the 3d level is filled do pairs with opposing spins form. Therefore. each 
atom in a transition metal has a permanent magnetic moment. dependent in strength to 
the number of unpaired electrons. Each atom behaves as a magnetic dipole. [Ref. 14: 
p. 609] 


B. DOMAINS 

Ferromagnetic behavior is due to unfilled energy levels in the 3d level (for iron, 
nickel, cobalt). Ferromagnetic materials have their powerful influence on magnetization 
because of the positive interaction between the dipoles of neighboring atoms. Within 
the grain structure of a ferromagnetic material, a substructure composed of magnetic 
domains is produced, even in the absence of an external field. Domains are regions 
which all of the dipoles are aligned to a common direction. In a material that has not 
been exposed to an external magnetic field, the individual domains have a random ori- 


entation and the net magnetization 1s zero. [Ref. 14: p. 613] 
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Fisure 34. Electron Spins of Transition Elements. 


The concept of a domain is key to understanding magnetism in materials. Domains 
are typically very small, measuring about .005 cm or less. Bloch walls, which are 
boundaries that separate the individual domains measure about 1000 A thick. These 
Bloch walls are narrow zones in which the direction of the magnetic moment gradually 
and continuously changes from that of one domain to that of the next (Figure 35 on 
pace G7). [her Iseor on, 

When an external field is anphed to tne material, domains that are alreadv lined up 
with the field grow at the expense of unaligned domains. In order for the domains to 
grow, the Bloch walls must move. The imposed magnetic field provides the force re- 
quired for the walls to migrate. As the strength of the field increases, favorably onented 
domains continue to grow and a greater net magnetization occurs. The saturation 
magnetizauon, produced when all of the domains are properly oriented, 1s the greatest 


amount of magnetization that the material can obtain. [Ref. 14: p. o14] 


C. HYSTERESIS LOOP 

When the field is removed, the resistance offered by the domain walls prevents re- 
growth of the domains into random orientations. As a result, many of the domains re- 
main oriented near the direction of the original field and a residual magnetization, 
known as the remanence B, persists in the material. The material now acts as a perma- 


nent magnet. [Ref. 14: p. 614] 
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Figure 35. Magnetic NIoments in Adjoining Atoms. 


If we now applwv a field in the reverse direction, the domains grow with an alignment 
in the opposite direction. A coercive force H, is required to force the domains to be 
randomly oriented. Further increases in the strength of the field eventually ahgn the 
domains to saturation in the opposite direction. [Ref. 14: p. 615] 

As the field continually alternates, the magnetization versus field relationship traces 
out a hysteresis loop (Figure 36 on page 68). This type of curve applies to all magnetic 
Mietenialss) alenals “hich have a low coercive force are low energy or soft magnetic 
materials. High coercive forces occur in rare earth permanent magnets. The power of 
the magnet is related to the size of the hysteresis loop, or the maximum product of B 
and H along the B-H curve. The area of the largest rectangle that can be drawn in the 

econd or fourth quadrants of the B-H curve is related to the energy required to 
demagnetize the magnet. Therefore, large hysteresis loops are desired for permanent 
magnets. as the loop also represents the amount of energy the magnet can store. [Ref. 
14: p. 617] 


D. MAGNET DEVELOPMENT 
Development of strong permanent magnets, often said to be magnetically hard, 1s 
aimed at improving both the magnetic permeability and the stability of the domains. 


Improved permanent magnets are produced by making the grain size so small that only 
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Figure 36. The Ferromagnetic Hysteresis Loop: showing the effect of the mag- 
netic fieid on inductance or magnetization. The dipole alignment leads 
tO SAlUration masnetizauon. poMmt ss, a Tenianence, point 4) audeo 
Coercive force, pene: 


one domain 1S present in each grain. Now the boundaries between domains are grain 
boundaries rather than Bloch walls; the domains can only change their orientation by 
rotatine. Which requires greater energy than domain growth. {Ret 14: 5. Geil] 

Two techniques are used to produce these magnetic materials, phase transf- 
ormations and powder metallurgy. Alnico, one of the most common of the complex 
metallic allovs, has a single-phase BCC structure at high temperatures. But when Alnico 
Slowly cools below $00 °C, a second BCC phase rich in iron and cobalt precipitates. 
This second phase is so fine that each precipitate particle is a single domain, producing 
a very high remanence, coercive field, and power. Often the alloys are permitted to cool 
and transform when in a magnetic field to align the domains as they form. [Ref. 14: p. 
621] 
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A second technique. powder metallurgy is used for another group of metallic allovs, 
including samarium-cobalt. The brittle intermetallic is crushed and ground to produce 
a fine powder in which each particle is a domain. The powder is then compacted when 
in an imposed magnetic field to align the powder domains. Careful sintering to avoid 


growth of the particles produces a solid powder metallurgy magnet. [Ref. 14: p. 621] 


B BENMPERATURE EFFECTS 

When the temperature of a ferromagnetic material is increased, the added thermal 
energy reduces the magnetic permeability or magnetization and permits the domains to 
become randomly oriented. Consequently, magnetization, remanence and the coercive 
field are all reduced at high temperatures. If the temperature exceeds the critical Curie 
temperature, magnetic properties of the material no longer exists. The Curie temper- 
ature corresponds to a threshold required to destrov the ordering of electron spins, so 
that the the permanent dipoles in each cell no longer exist. The Curie temperature de- 
pends on the tvpe of magnetic material and can be changed by alloying elements. Two 
common ferromagnetic materials are iron and cobalt. Each has a Curie temperature of 
tre C and 1131 °C respectively. [Ref. 14: p. 619} 


B. RADIATION EFFECTS 

Radiation effects have caused increasing concern in recent years. Experiments con- 
ducted on present commerciallv-available permanent magnets have shown no harmful 
effects on magnetic stability when magnets were subjected to moderate levels of radi- 
ation. Magnet technology is at its optimum, and with the improvement to produce very 
high remanence and coercive field, more and more applications of permanent magnets 


are in high radiation environments. 
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APPENDIX B. GLOSSARY OF TERMS 


The following list of common terms associated with magnetic properties of materials 
are provided to aid the reader. The publication Permanent Magnet Guidelines by Mag- 
netic Materials Producers Association is the source from which the following terms were 


Originated: 


B, magnetic induction, is the magnetic field induced by a field strength, H, at a given 
point. It is the vector sum, at each point within the substance, of the 
magnetic field strength and resultant intrinsic induction. Magnetic in- 
duction 1s the flux per unit area normal to the direction of the magnetic 


path. 


B, , remanent induction, is any magnetic induction that remains in a magnetic material 
after removal of an applied saturating magnetic field, H,. (B, 1s the mag- 
netic indcution at any point on the demagnetization curve; measured in 


gauss.) 


B,/H, , energy product, indicates the energy that a magnet material can supply to an 
external a magnetic circuit when operating at a point on its demagnetization 


curve; measured in mega-gauss-oersteds. 


(By/Ha)max » aAXimum energy product, indicates the maximum intrinsic induction pos- 


sible in a material. 


B, , intrinsic induction, 1s the contribution of the magnetic material to the total mag- 
netic induction, B. It is the vector difference between the magnetic in- 
duction in the material and the magnetic induction that would exist in a 
vacuum under the same filed strength, H. This relation is expressed by the 


equation: 


where B,; = intrinsic induction in gauss; B = magnetic induction in gauss; 


}Hi = field strength in oersteds. 
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B,, » remanent induction, 1s the magnetic induction that remains in a magnetic material 


after magnetizing and conditioning for final use; measured in gauss. 


B,,H, » energy product, indicates the energy that a magnet material can supply to an 
external magnetic circuit When operating at the point (B,,H,,); measured in 


mega-gauss-oersteds. 


B, , residual induction (or flux density), is the magnetic induction corresponding to 
zero magnetizing force in a magnetic material after saturation in a closed 


circult; measured in gauss. 


H, magnetic field strength, (magnetizing or demagnetizing force), is the measure of the 
vector magnetic quantity that determines the ability of an electric current, 
or a magnetic body, to induce a magnetic field at a given point; measured 


in oersteds. 


H, , coercive force, of a material is the demagnetizing force corresponding to zero 
magnetic induction, B, in a magnetic material after saturation; measured in 


oersteds. 


H,, , intrinsic coercive force, of a material, indicates its resistance to demagnetization. 
It is the demagnetizing force corresponding to Zero intrinsic induction in a 


material after saturation: measured in oersted. 


ally is that value of H corresponding to the remanent induction, B,; measured 


in oersteds. 


He , is that value of H corresponding to the remanent induction, B,,; measured 


in oersteds. 


T, , Curie temperature, is the transition temperature above which a material loses its 


permanent magnet properties. 


u, permeability, is a general term used to express various relationships between mag- 


netic induction, B, and the field strength, H. 


@, magnetic flux, is a measurable concept that evolved in an attempt to describe the 


“flow” of a magnetic field. 


closed circuit condition exist when the external flux path of a permanent magnet 1s 


confined within high permeability material. 


71 


demagnetization curve is the second (or fourth) quadrant of a major hysteresis loop. 


Points on this curve are designated by the coodinates B, and H,. 


fluxmeter is a galvanometer that measures the change of flux linkage with a search 


Colle 


hysteresis loop is a closed curve obtained for a material by plotting (usually to rectan- 
gular coordinates) corresponding values of magnetic induction, B, for 
ordinates and magnetizing force, H, for abscissa when the material is pass- 
ing through a complete cycle between definite limits of either magnetizing 


force, H, or magnetic induction, B. 


major hysteresis loop of a material is the closed loop obtained when the material is 


cvcled between positive and negative saturation. 


open circuit condition exist When a magnetized magnet 1s by itself with no external flux 


path of high permeability material. 


permanent magnet is a body that is capable of maintaining a magnetic field at other 


than cryogenic temperatures with no expenditure of power. 


magnetic saturation of a material exists When an increase in magnetizing force, H, does 


not cause an increase in the intrinsic magnetic induction, B, of the material. 


search coil is a coiled conductor, usually of known area and number of turns, that is 


used with a fluxmeter to measure the change of flux linkage with the coil. 
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APPENDIX C. ITS OVERVIEW 


The Integrated Tiger Series of Electron’Photon Transport codes (ITS) is the most widely 
used particle transport code in the world. The code package was developed to incorpo- 
rate eight individual codes which were developed over the period from 1968 to 1981. 
All the codes are based on the original ETRAN model developed by M. Berger and D. 
Seltzer. The ITS code system consists of four primary code packages [Ref. 12: p. 9]: 


XMDATA: The electron and photon cross section data file. 
ANGEN: The cross section generation program. 
ITS: The Monte Carlo program file. 
UCPEML: A machine portable update emulator. 


The heart of the ITS is the program library file ITS, which contains the eight Monte 
Carlo programs plus system directives for the CRAY, IBM, VAX, and CDC operating 
systems. The update emulator program UPEML creates the various Monte Carlo codes 
for a given svstem with any corrections to those codes that may be desired. The output 
fortran source code from LPEML is then compiled and stored as an executable module. 
Program NGEN generates the problem specific cross section data tape using file 
NDATA for referenced inputs and a user defined input file. The Monte Carlo codes 
then read in the cross section tape and process the user defined problem. 

One of the eight ITS codes is CYLTRAN, which simulates the transport of particle 
trajectories through a three-dimensional multimaterial cylinder. For this project only the 
CYLTRAN code was required. As an ITS user the following steps are required to exe- 


cute an ITS run: 


1. Generate a cross section tape based on the different type of materials contained in 


the cylindrical geometry of a problem. 


2. Create an input file which list all the input parameters required to calculate desired 


outputs. 


ted 


Submit the input file and the generated cross section tape to the ITS Monte Carlo 


codes to execute a run. 
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Table 12 on page 74 1s a sample input file to generate the cross section tape for the 
materials in a cvlindical geometry. Each material line represents a different medium in 
the cylinder. Percentages of each material in a compound and its density must be spec- 
ified. Single element lines such as Cu, has its density stored in ITS and 1s automatically 


entered into the simulation. 


Table 12. SAMPLE INPUT FILE TO CREATE A CROSS SECTION TAPE 


Energy 79.0 
Material Cu .04 Ni .07 W .89 Density 17.1 
Material Cu 
Material N .80 O .20 Gas Density .001226 


Material Li .72 F .28 Denstiy 2.69 
EME 
79 NleV Cross Section for Prod Tet, Air and LiF 





Once a cross section tape is generated, an input file with the parameters design for 
a particular simulation must be created. Table 13 on page 76 is a sample input file to 
execute an [TS run. The keywords at the beginning of each line are relatively self- 
explanatory. However, few kevwords and their numerical parameters require some ex- 
planation to fully understand their importance. These few keywords and their 


explanation are listed [Ref. 12: pp. 21-23]: 


Geometry This keyword sets up the cylindrical geometry of the problem into input 
zones, Where in this case there are seven input zones. Each line of numerical 
parameters following the keyword describes the dimensions of each zone. 
Figure 37 on page 77 1s a graphic representation of the seven input zones. 
Note that input zone seven is divided into 10 subzones. This figure repres- 
ents the half section of a problem cylinder where the z-axis is the axis of the 
beam center. Figure 37 on page 77 is also the geometrical model used at 
LLNL. The fifth column of numerical parameters corresponds to the ma- 


terial line as listed in the sample input file in Table 12. 
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Electron/Photon-Escape This keyword tallies the number of incident electrons and 
photons that escapes the cylindrical geometry, either latterally or transmit- 


ted through. 
NBINE _ tallies the escaped electron’photons in specified energy bins. 
NBINT _ tallies them in angular bins. 


Electron/Photon-Flux This keyword tallies the energy deposition of electrons and 


photons in the subzones described in the kevword geometry. 


Histories This keyword represents the number of primary particle histories to follow 


as it transport through each medium in the cylinder. 


Using the input file in Table 12 to generate a cross section tape and executing the 
input file in Table 13 on page 76 with this tape, a simulation run measuring gamma and 
electron dose deposited on L1,F material (TLD chips used at LLNL) will be executed. 
Table 14 on page 78 1s a sample output file used to plot radiation dosage patterns of 
various materials. This file was generated from the energy deposition output section 
mreated by ITS. 


Table 13. SAMPLE OF AN INPUT FILE (0 EGE Its 


Echo | 
TITLE 
Mat:Tet,Cu,air, LIF R=.75cm D=78.7cm His:10000 E=79MeV 
Electrons 
Energy 79.0 
Cutoffs .1 .1 
Position 0.0 0.0 0.0 
Radius 0.75 
Direction 0.0 0.0 
Geometry 7 


0.000 0.254 0.000 0.944 I 
0.000 0.254 0.944 3.492 2 
0.009 0.254 3.492 10.000 3 
0.254 1.054 0.000 S.B2 = 
0.254 1.054 3.492 10.000 R 
1.054 78.709 0.000 10.000 3 
78.700 79.060 0.000 10.000 4 1 1 = 10 


Electron-Escape 
NBINE 23 User 
75. 70. 60. 55. 50. 45. 40. 35 
25. 205 15. 1058.7 eoeo. 4.3. 
6.) 
NBINT 16 User 
1.0 2.0 3.0 4. 5. 6. 7. 8. 9. 10. 
14. 20. 30. 45. 90. 180. 
Photon-Escape 
NBINE 23 User 
75. 70. 60. 55. 
25.20, Isai: 
6.1] 
NBINT 16 User 
1.0 2.0 3.0 4. 5.6. 72.879 
15. 29. 30. 45. 90. 180. 
Electron-fiux 1 10 
NBINE 10 
Photon-Flux | 10 
NBINE 10 
Histories 10000 


: we) 
WZ 


0. 4 
ql. 


ry te 
_ 
we) 
re) 
Ts 


rags 
typ 


co 
nu 


76 


79.9 


Coa 





Distance (cm) 


<= 
Lo) 
= 
+ 
Ve) 
~ 
© 
= nN 
S $3 . 
(wo) snipey | 
Figure 37. Graphic Representation of the Keyword Geometry: where 


heavimet= 1, copper=4, air=3,5,6 and the absorbing material= 7. 
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Table 14. SAMPLE OF A RADIATION DOSAGE PATTERN OUTPUT FILE 


DISTANCE PRIMARIES |= GAMMAS TOTAL 
(CM) (GRads/Coul) | (GRads/Coul)  (GRads/Coul) 
0.5000 0.0114 54 0.0648 25 0.0835 20 
1.5000 0.0081 56 0.0470 14 0.0557 15 
2.5000 0.0165 22 0.0377 11 0.0535 12 
3.5000 0.0244 13 0.0430 12 0.0657 7 
4.5000 0.0260 16 0.0269 12 0.0527 9 
5.5000 0.0183 29 0.0238 8 0.0419 10 
6.5000 0.0161 10 0.0227 12 0.0374 9 
7.5000 0.0169 20 0.0211 9 0.0381 8 
8.5060 0.0194 9 0.0188 9 0.0371 4 
9.5000 0.0185 7 0.0150 7 0.0315 6 
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